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A  Mathematical  Model  for  the  Class  V 


Abstract 


The  purpose  of  this  Investigation  is  to  develop  a  mathematical 
model  for  the  Class  V  f iextensional  underwater  acoustic  transducer. 

The  transducer  is  approximated  through  the  consideration  of  three 
distinct  problems.  The  problem  of  a  thin  piezoelectric  disk  with  an 
arbitrary  impedance  on  its  edge  is  solved  in  terms  of  Bessel  functions. 
The  shell  vibration  problem  is  solved  using  a  finite  difference  model 
to  approximate  the  shell.  The  acoustic  radiation  problem  is  solved  by 
obtaining  the  source  density  distribution  for  \  system  of  quadri¬ 
laterals  representing  the  transducer.  With  the  source  density  of  each 
quadrilateral,  the  near  and  far  field  pressures  and  velocities  can  be 
found.  Utilizing  these  three  components,  a  model  is  then  constructed 
for  the  transducer. 

A  comparison  of  the  results  from  the  mathematical  model  with 
those  obtained  from  experiments  is  made  in  order  to  validate  the  model. 
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1.  INTRODUCTION 


In  general,  the  various  f lextenslonal  designs  can  be  placed  in 
one  of  five  different  classes  [2].  The  purpose  of  this  investigation 
is  to  develop  a  mathenatical  model  for  the  Class  V  f lextenslonal 
underwater  acoustic  transducer.  A  picture  of  a  Class  V  transducer  is 
shown  in  Figure  1.1  and  a  typical  sketch  in  Figure  1.2.  The  Class  V 
flextensional  transducer  design  was  originally  proposed  as  having  pos¬ 
sible  applications  as  a  sonobuoy  transducer. 

As  shown  in  Figure  1.2,  a  Class  V  design  consists  of  two  shallow 
spherical  shells  bonded  at  the  common  boundary  by  an  epoxy  cement  to 
the  edge  of  a  chin  piezoelectric  disk.  The  piezoelectric  disk  is  iso¬ 
lated  electrically  from  Che  two  shells  by  the  removal  of  the  electrodes 
from  the  region  where  there  is  contact  between  Che  disk  and  Che  shells. 
Sufficient  epoxy  is  applied  so  as  to  firmly  attach  each  of  the  shells 
CO  Che  disk.  Two  small  holes  are  drilled  through  Che  shells  and  serve 
as  entrances  for  the  electrical  leads  of  Che  disk. 

To  attain  Che  stated  objective,  one  develops  a  mathematical  model 
for  each  of  the  important  components  of  the  system.  Then,  combining 
them  in  a  manner  similar  to  that  used  by  Royster  [21],  one  obtains  an 
approximate  working  model  for  the  coi^>leCe  system.  The  model  is  then 
programmied  on  the  IB1  360-75  computer,  and  an  evaioi.(ion  is  made  by 
comparing  experimental^  and  predicted  results. 

^Dahlke,  H.  E.  and  L.  H.  Royster.  1964.  Unpublished  notes  on 
Shallow  Shell  Transducer.  North  American  Aviation,  Inc..  Columbus, 
Ohio. 


4 


2.  REVIEW  OF  LITERATURE 

2. 1  General 

The  developmeni:  of  detailed,  complete  mathematical  models  for 
the  various  classes  of  f lextensional  underwater  acoustic  transducers 
has  been  underway  for  about  five  years.  Hence,  the  literature  in 
this  area  is  limited  to  only  a  few  papers  and  technical  reports.  On 
the  other  hand,  elements  of  the  transducer  problem,  such  as  the  thin 
piezoelectric  disks,  spherical  caps,  and  acoustic  radiation  problems 
for  different  geometries  have  been  studied  extensively. 

2.2  Mathematical  Models  of  Flextensional  Underwater 
Acoustic  Transducers 

The  present  state  of  the  art  of  analytic  models  for  flextensional 
underwater  acoustic  transducers  is  presented  by  Brigham  and  Royster  [2]. 
They  also  classified  the  various  flextensional  transducer  designs  into 
five  different  classes.  Presently,  there  are  two  analytic  models,  one 
complete  and  one  partial,  which  have  been  presented.  The  complete 
model  was  presented  by  Royster  [21]  for  the  Class  I  transducer  as 
shown  in  Figure  2.1.  The  partial  model  was  presented  by  Boone  and 
Royster  [1]  for  the  Class  II  transducer  as  shown  in  Figure  2.2.  Both 
models  use  a  finite  difference  model  for  the  shell  but  differ  in  the 
manner  by  which  the  shell  and  stack  are  coupled  and  in  the  media  in 
which  they  operate.  The  Class  I  model  assumes  that  the  shell  vibrates 
Tndependently  of  the  stack  and  then  solves  the  eigenvalue  problem  for 
the  shell.  On  the  other  hand,  the  Class  II  model  couples  the  electrode 
shorted  stack  to  the  shell  through  the  edge  segment  equation  of  motion 


7 


-I 

1 


and  then  solves  the  eigenvalue  problem.  The  Class  I  model  can  also 
operate  in  media  which  produce  acoustic  loads  while  the  Class  II  model 
in  its  present  state  of  development  is  limited  to  media  which  produce 
only  negligible  acoustic  loads. 

Presently,  work  is  underway  on  discrete  models  for  Classes  II, 

III  and  IV  at  the  Center  for  Acoustical  Studies,  and  a  continuous 
closed  form  solution  for  a  Class  IV  oval  is  being  developed  by  North 
American  Aviation. 

Since  the  mathematical  model  to  be  developed  herein  is  for  the 
Class  V  transducer,  a  review  of  the  literature  in  the  areas  of  thin 
piezoelectric  disks,  spherical  caps,  and  acoustic  radiation  problems 
for  three-dimensional  bodies  is  of  interest. 

2.3  Thin  Piezoelectric  Disks 

Three  works  in  the  literature  are  pertinent  to  an  investigation 
of  the  vibration  of  a  thin  piezoelectric  disk.  The  first  of  these  is 
the  vibration  of  a  free,  thin,  electrostrictive  disk  by  Mason  [11]. 
Mason  established  the  necessary  system  of  equations  and  the  assump¬ 
tions  and  procedure  needed  to  solve  the  problem.  The  resulting  dis¬ 
placements  differed  from  the  elastic  solution  for  a  thin  disk  in  that 
the  material  properties  were  measured  in  a  constant  electric  field. 

The  second  work  by  Tachibana  [23],  used  Mason's  procedure  and 


solved  the  problem  of  a  partially.plated,  thin  piezoelectric  disk. 
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frequencies,  and  dynamic  electromechanical  coupling  coefficient  for  a 
thin  piezoelectric  disk  with  an  arbitrary  impedance  on  its  boundary. 
The  expressions  were  specialized  for  the  case  of  a  spring  like 
boundary  condition. 

2.4  Spherical  Caps 

Langhaar  [10]  developed  an  expression  for  the  strain  energy  in 
thin  elastic  shells.  Special  cases  for  the  geometry  were  then  c(n- 
sidered,  one  of  the  cases  being  the  sphere.  Froff  Langharrs  strain 
energy  expressions  for  a  sphere,  McDonald  [12]  obtained  the  equations 
of  motion  for  a  spherical  cap  with  a  clamped  edge,  using  finite  dif¬ 
ferences.  These  equations  were  then  solved  numerically  for  the 
eigenvalues  and  eigenvectors  of  the  free  undamped  problem  and  were 
compared  to  other  analytic  results.  The  forced  vibration  problem  was 
also  solved  for  various  loads. 

Using  the  technique  presented  in  his  1959  work,  McDonald  [13] 

Chen  presented  the  equations  of  motion  for  free  undamped  vibration  of 
a  thin  shell  of  revolution  composed  of  an  orthotropic  material  with 
variable  elastic  properties  and  chickness. 

Nelson  and  Royster  [15],  utilizing  the  techniques  presented  in 
both  papers  by  McDonald,  derived  the  equations  of  motion  for  the 
guided-pinned  and  guided-clamped  spherical  cap.  The  results  were  then 
set  up  in  matrix  representation  in  order  to  obtain  a  numerical  solu- 
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2.3  Acoustic  Radiation 

Three  numerical  techniques  were  considered  as  possible  methods 
for  obtaining  the  radiation  field  for  the  Class  V  flextensional 
transducer.  This  section  will  consist  of  a  review  of  those  techniques. 

Chertock  [3]  developed  a  method  for  determining  the  radiation 
field  of  a  vibrating  body.  In  general^  a  Fredholm  integral  equation 
is  solved  numerically  to  give  surface  pressures.  Then  the  special 
case  of  a  surface  of  revolution  with  an  arbitrary  velocity  distribu¬ 
tion  and  frequency  is  described  in  detail. 

Hess  [8]  developed  a  method  of  solving  the  radiation  problem  for 
an  arbitrary  three-dimensional  body  represented  as  a  system  of 
quadrilaterals.  Specification  of  the  geometry^  normal  surface 
velocities,  and  wave  number  leads  to  a  set  of  linear  algebraic  equa¬ 
tions  which  approximate  an  integral  equation  for  the  source  density 
distribution.  Solution  of  the  system  gives  near  and  far  field  pres¬ 
sures  as  well  as  other  information. 

Schench  [22]  developed  four  different  integral  formulations  for 
obtaining  approximate  solutions  to  the  exterior  steady-state  acoustic 
radiation  problem  for  an  arbitrary  surface  whose  normal  velocity  is 
specified.  The  fourth  method,  a  Combined  Helmholtz  Integral  Equation 
Formulation  (CHIEF),  overcomes  numerical  difficulties  that  arise  for 
certain  wave  numbers. 

The  method  chosen  for  use  herein  is  that  developed  by  Hess.  This 
method  was  chosen  due  to  its  general  approach  to  the  body  shape.  The 
limitations  of  the  method  in  the  frequency  range  presented  no  problem 
since  the  fundamental  mode  of  the  Class  V  was  well  below  the  first 


forbidden  frequency." 
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3.  DEVELOPMENT  OF  MATHBIATICAL  MODEL 

3. 1  General 

To  develop  a  mathematical  model  for  the  Class  V  flextenslonal 
transducer^  the  model  Is  divided  into  three  sections.  The  first  part 
is  the  vibration  of  a  thin  piezoelectric  disk  with  an  arbitrary 
impedance  on  its  edge  [14].  The  second  is  the  vibration  of  a 
spherical  cap  having  a  horizontally  guided-pinned  boundary  condition 
[15].  The  third  section  is  the  acoustic  radiation  problem  due  to  the 
surface  displacements  on  the  shell. 

The  manner  in  which  these  three  problems  are  joined  depends  on 
the  assumptions  made.  Hence,  discussion  of  this  aspect  of  the  overall 
transducer  model  will  be  delayed  until  after  the  three  problems  have 
been  considered. 


3.2  Equations  of  State 

Both  the  equations  of  motion  and  the  linear  adiabatic  equations 
of  state  for  a  piezoelectric  material  in  a  rectangular  cartesian 
coordinate  system  were  outlined  by  Royster  [20]  and  the  practical 
engineering  limitations  noted.  Therefore,  only  the  resulting  equa¬ 
tions,  along  v.'ith  the  assumptions  made,  will  be  noted  here.  In 
general,  the  notation  used  will  also  follow  that  outlined  in  the 
IRE  (1949)  publication  [17]. 

The  equations  of  motion  arise  out  of  the  application  of  Newton's 


d^u. 


O.l) 


dt 


where  T.  is  a  symmetric  tensor  %rhen  body  moments  are  neglected. 

^  J 
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T 


On  the  other  hand^  the  equations  of  state  r:annot  be  so  easily 
determined.  First  using  the  conservation  of  energy  and  then  assuming 


that  the  dissipation  of  energ5r^by  the  material  is  small  compared  to 
the  rate  of  change  of  kinetic  and  potential  energy,  and  that  all 
processes  are  reversible,  one  can  obtain 


\  -  "  Vl  *  V 


(3,2) 


where  (  )  denotes  differentiation  with  respect  to  time.  From  equation 
(3.2),  it  can  then  be  assumed  that 


%  =  “i'  • 


<3.  3) 


From  equation  (3.3),  one  can  derive  the  two  equations  of  state  needed 
in  the  development  of  the  mathematical  model.  In  obtaining  these 
equations  of  state,  one  must  assume  the  following; 

1)  small  deflection  gradients  (linear  theory), 


2)  adiabatic  process  (no  heat  transferred). 

These  assumptions  limit  the  transducer  to  low  power  applicatinus .  The 
resulting  equations  of  state,  which  are  cooinonly  called  the  lintic 
adiabatic  equations  nf  state,  are  given  in  rectangular  cartesian  tcn*^or 
form  by 


(3.4) 


+  e  E  . 
nn  B 


D  =  d  ,  T. 
n  nki  ki 


(3.5) 
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Equations  (3.4)  and  (3.5)  will  be  referred  to  simply  as  equations  of 
state. 

As  shown  in  Figure  3.1,  the  coordinate  system  used  for  the  thin 
piezoelectric  disk  is  a  cylindrical  coordinate  system.  Use  of  the 
appropriate  tensor  transformations  and  material  properties  as  shown 
in  Appendix  7. 1  results  in  the  equations  of  state  for  cylindrical 
coordinates  given  by 


S  =  S^.T  +  +  T  )  +  d,-E 

rr  11  rr  12  98  **  31  z 


S  +  S®  (T  +  T  )  +  d^.E 

09  11  68  12  rr  zz  31  z 


S  =  sf.T  +  s!‘-(T  +  T 

zz  11  zz  12  rr 


^re  ■  ^^11  ■  ^12^'^re 


S  =  (sf  -  )T 
rz  11  12  rz 
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3.3.1  Electric  Field  and  Piezoelectric  Equations 

The  disk  being  considered  is  a  thin  circular  disk  with  electrodes 
on  the  top  and  bottom  as  shown  In  Figure  3.2.  The  thickness  of  the 
dlsk^  h^^  Is  much  less  than  the  radius  of  the  disk,  a^, 

"d « ‘d- 

In  general,  the  displacements  In  cylindrical  coordinates  are 
denoted  by 

Uj.  = 

Ug  =  u^(r,  0,  z,  t)  (3.  7) 

and 

“z  ""  • 

Now  the  motio|^will  be  assumed  to  be  primarily  radial.  Therefore, 
there  will  be  no  displacement  allowed  in  the  0  direction,  and  all 
other  displacements  must  be  Independent  of  0.  Also,  the  displacement 
in  the  radial  direction  will  be  required  to  be  independent  of  z. 
Equation  (3. 7)  can  now  be  written  as 

“z  " 

and 


The  stresses  on  a  cylindrical  clement  are  shown  in  Figure  3.3. 
The  bottom  and  top  of  the  disk  are  free  surfaces.  Hence,  since  the 
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disk  is  thin^  one  can  assume  that 


T  =  T.  =  T  =  0  . 
rz  6z  zz 


(3.9) 


The  strain-displacement  relations  for  cylindrical  coordinates  are 
given  by  Fung  [5]  to  be 


rr 


Su 

r 

w 


66  r 


1 


ST 


zz 


Su 

z 

(5z  ^ 


2S 


6  ^ 


u 


«  1 

6,1  r 

’ 


(3.10) 


2S 


rz 


^u  Su 

W-  'ST 


2S 


,  du 
1  z 


du. 


6z 


Substitution  of  equation  (3.8)  in  equations  (3.10)  yields 


rr 


2S 


rz 


Su 

r 


Su 

z 

dr 


^66  r  ^ 


zz 


au 

z 

^z  ’ 


and 


^6  =  'e^  = 


(3.11) 


The  equilibrium  equations  in  cylindrical  coordinates,  which  can 

I 

be  obtained  from  equation  (3.1)  through  the  procedure  outlined  in 
Appendix  7.1,  are  given  by 


Bt  ,  Bt  ,  Bt  t  -t„^ 

_  rr  ^  1  re  ^  rz  ^  rr  66 

-  dr- ^7  — F — 


„  1  aT„„  BT„  T  „ 

^  +  ez  .  „  re 

■■  dr~  ^  7  dr~  ^  ^  T" 


Bt  ,  Bt^  Bt  t 

_  rz  1  02  .  zz  rz 

p^'z "  dr~  7  dr~  dr~  ^  -r 


(3.12) 
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and 
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From  equations  (3.11)  and  (3.17),  one  can  note  that 


S 

rz 


=  0 


so  that 


(3.18) 


(3.19) 


and 

T  ^  =  0  .  (3.20) 

r6 

One  can  now  note,  from  equation  (3.20)  and  the  independence  of 
the  problem  with  respect  to  9,  that  the  equilibrium  equation  given  by 
(3.14)  is  satisfied  identically.  From  equations  (3.8),  (3.15),  and 
(3.19),  one  observes  that  u^  is  now  restricted  to  being  a  function  of 
z  and  a  linear  function  of  time.  Therefore,  it  is  now  assumed  that  u 
is  zero,  so  that  equation  (3.8)  becomes 

Uj.  =  u^(r,t)  (3.21) 


and 


u„  =  u  =0 

0  z 


The  strain-displacement  equations  given  by  (3.11)  are  reduced  to 


rr 


3u 

r 


and 


ee 


r 

r 


(3.22) 


The  equilibrium  equations  become 


St  t  -t^„ 

rr  .  rr  99 
^  ^  - ? - 


(3.23) 
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From  equations  (3.17)^  the  equations  of  state  for  the  two-dimrnsional 
radial  motion  of  the  disk  are  given  by 


S  =  S?,T  +  +  d-,E 

rr  11  rr  12  06  31  z 


(3.24) 


and 


(3.25) 


‘'3l(\r  "«33"z 


(3.26) 


3.3.2  Equation  of  Motion 

In  order  to  obtain  the  equation  of  motion,  one  must  find,  from 
equations  (3.24)  and  (3.25),  the  stresses  in  terms  of  strain  and 
electric  field.  From  Mason  [11],  it  is  round  that 


=  Y 


and 


’12 


=  a 


11 


11 


(3.27) 


Hence,  equations  (3.24)  and  (3.25)  yield 


Y^ 

T  =  (S  +  aS)  .  d..  ~E 

rr  ,  2  rr  00  31  1-a  z 

i-CT 


(3.28) 


ana 


Y^  Y^ 

(S„„  +  aS  )  -  d,,  E  . 
90  ,  2  '  60  rr'  31  l-o  z 

i-a 


(3.29) 


Substitution  of  equations  (3.28)  and  (3.29)  into  equation  (3.23)  gives 


2 

Y^  s  u  au^  u 

"  7T  F  a? - 2^ 

1-0  dr  r 


(3.30) 


Equation  (3.30)  is  the  equation  of  motion  for  a  thin  piezoelectric 
disk.  This  is  identical  to  the  equation  of  motion  for  a  thin  elastic 
disk  and  electrostrictive  disk  [11],  except  for  the  fact  that  Young's 
modulus  in  this  case  is  measured  under  a  constant  electric  field  con¬ 
dition.  Equation  (3.30)  is  also  identical  to  that  obtained  by 
Tachibana  [23],  as  it  should  be. 


3.3.3  Boundary  Conditions  and  Initial  Conditions 

The  model  developed  herein  is  limited  to  steady  state  vibrations; 
therefore,  initial  conditions  are  not  necessary.  Thus,  one  assumes  a 
steady  state  displacement  represented  by 

u^(r,t)  =  U^(r)ej'"*^  .  (3.31) 

The  first  boundary  condition  to  be  given  for  the  disk  applies  to 
its  center.  Since  the  only  motion  allowed  is  radial,  then  the  center 
must  not  move;  therefore, 

Uj.1  =0.  (3.32) 

r=o 


On  the  edge  of  the  disk,  an  arbitrary  specific  impedance  is  applied. 
Hence,  from  Figure  3.4  and  the  definition  of  specific  impedance  [91, 
one  obtains 


n  +  jr  . 


(3.33) 


The  arbitrary  impedance  consists  of  a  real  (resistive)  part  and  an 
imaginary  (reactive)  part. 
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To  obtain  the  displacements  of  the  disk,  one  first  ases  the 
steady  state  assumption  given  by  equation  (3.31)  in  the  equation  of 
motion  given  by  equation  (3.30)  to  obtain 


dU 

r 

sr 


+  (x^-l)U^ 


=  0  , 


(3.34) 


where 


and 


(3.35) 


=  Y^/  p(l-a^)  . 


(3.36) 


The  parameter  x  will  be  called  the  dimensionless  frequency  parameter, 
while  V  is  the  wave  speed  in  the  disk. 

The  solution  of  equation  (3.34),  which  is  a  form  of  Bessel's 
equation,  is  given  in  Churchill  [4].  Hence,  one  can  write 

-  AJ^(x)  +  BY^(x)  ,  (3.37) 

where  A  and  B  are  complex  constants. 

The  boundary  conditions  given  by  equations  (3.32)  and  (3.33)  are 
now  used.  In  order  to  satisfy  (3.32),  one  must  require  B  -  0  so  that 
equation  (3.37)  becomes 

=  AJj(x)  .  (3. 38) 

Stresses,  strains  and  the  electrical  field  are  also  steady  state  forms. 


i 
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For  purposes  of  convenience,  their  exponentials  will  be  dropped. 
Therefore,  expressions  and  equations  will  be  amplitude  relations  un¬ 
less  otherwise  stated.  Using  equation  (3.28)  in  (3.31)  and  the  result 
from  this  substitution  in  equations  (3.22),  the  steady  state  strains 
are  found.  The  strains  may  then  be  used  in  equations  (3.28)  and 
(3.29)  to  obtain  the  steady  state  stresses.  In  particular,  the 
radial  stress  is  found  to  be 


r  .A 

rr  ,  2 

1-0 


.E  V. 


J 1  (^)  V"  '  n 

.  (i.a) 

a 


(3.39) 


where  a  steady  state  voltage  of  the  form 


V  =  V  e 
o 


j(i)t 


(3.40) 


has  been  assumed.  Hence,  using  equations  (3,31),  (3.38)  and  (3.39)  in 
the  boundary  condition  given  by  equation  (3.33),  one  can  obtain  the 
complex  constant  A  so  that  the  radial  displacement,  equation  (3.38), 
becomes 


U^(r) 


Y(TVj6)  -^1^*^ 


T]^  +  8^ 


(3.41) 


where 


31  l-o  h. 


(3.42) 


n 


(3.4:) 
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,  _  ,„2  ’'Si  v(T(.i6)  ■'l^V 
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2  F 

2d  V 

2,  T  ,  o 

^  ^  ~T~r^  HT  • 

d 


(3.55) 


Using  equation  (3.55)  in  (3.52),  the  electrical  admittance  is 


given  by 


Y  =  g  +  jb  , 


(3.56) 


where 


g  -  2™a2  ^^"^31  V6 

®  -  ™^d  ^2  ^  ^2 


=  2  r  ^^31  Vn 

™  d  V^(Ua)v 


(3.57) 


Si 
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(3.58) 


which  is  the  planar  coupling  coefficient  as  given  in  the  IRE  standards 
of  1958  [18].  From  equations  (3.51)  and  (3.56),  the  electrical 
impedance  can  be  shown  as 


7  =  8  -  jb 

2^,2 

g  +b 


(3.59) 
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3.4  Shell  Vibrations 

As  stated  previously,  the  vibration  of  a  spherical  cap  having  a 
horizontally  guided-pinned  boundary  condition  has  been  studied  by 
Nelson  and  Royster  [15].  Their  work  was  presented  in  enough  detail 
so  that  a  basic  review  will  be  all  that  is  presented  here. 

3.4.1  Strain  Energy  in  a  Sphere 

From  the  results  presented  by  Langhaar  [10],  the  strain  energy 
due  to  stretching  and  herding  can  be  obtained.  The  strain  energy  due 
to  stretching  is  given  by 


2tt  e 


U.  =1-^  r  r  (  +CSC  e(v„  +  u  cos  6 +w  sin  6) 

i  i-O*  J  J  0  vp 


O  0 


+  2ff  CSC  8(u  +w)  (v  +  u  cos  0  +  w  sin  0) 
6  8 


+  csc^0(u  +  V  sin  0  -  V  cos  0)  ^  }  h  sin9d6dcp  (3.60) 

2  cp  0  s 


and  the  strain  energy  due  to  bending  by 


U 


2  12(l-cr) 


2tt  0  2  4 


o  o 


+  2CTCSC  0(Wgg)(WgSin  0cos  0 +W5p5p) 


+  2(l-o)csc^9(w  cot0-w.  )^)/  (h^/a^)  sin  0 d9  dtp 
Cp  9cp  s  s  ^ 


(3.61) 


where  ~  ^  ■  •  •  The  coordinates  system  used  for  the 

spherical  cap  is  .shown  in  Figure  3.5.  The  potential  energy  of  the 
external  forces,  can  be  readily  found  from  Figure  3.5  to  be 


given  as 


U 


ext 


-a^  J  J  (Xu  +  Yv  +  Zw)  sin  ededcp 


(3.62) 


o  o 


Thus,  the  total  strain  energy  V  of  the  system  is 


V  =  UjL  -H 


(3.63) 


Due  to  the  symmetrical  geometry  of  the  transducer,  one  is  only 
interested  in  rotationally  symmetric  vibrations.  Requiring  rotational 
symmetric  vibrations  and  thus  independence  with  respect  to  the 
circumferential  variable,  cp,  equations  (3.60),  (3,61)  and  (3.62) 
reduce  to 


Q  0  9  \  ‘ 

(  (liQ+w)^  9(u  cos' e +w  sin  0) 


2a CSC  e(u. +w)(ucos  e+wsin  6)  )  h^  sin  ©de  , 

U 


(3.64) 


^^2  = 


U  / 

^rru  ^  J  [  ,^^^  +  csc‘^0(w^sin  ecos  0)' 


12(l-a)ag  o 


2a  csc^0(WgQ)  (WgSin  ©cos  0)  }  hj  sin  0d6  , 


(3.65) 


and 


ext 


=  _2TTa^  r  (Xu  +  Zw)  sin  ed0 

S  V 


(3.66) 
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3,4.2  Boundary  Conditions 

Two  sets  of  boundary  conditions  are  needed  for  the  shell.  The 
first  set  describes  the  crown  point;  the  second  set  describes  the 
lower  edge  of  the  shell,  i.e.  9  = 

Rotationally  symmetrical  vibrations  about  the  axis  of  revolution 
have  been  required  for  the  shell.  Therefore,  the  only  possible  motion 
at  the  crown  point  is  one  in  which  the  crown  is  displaced  in  a  radial 
direction.  No  tangential  motion  can  accompany  the  radial  notion,  and 
from  geometric  considerations  the  slope  remains  zero.  These  require¬ 
ments  are  expressed  mathematically  as 


9=0 


=  0 


(3.67) 


and 


McDonald 


Sw 


>3 
o  w 


^9 


=  0  . 

9=0 

adds  to  this  set  the  requirement  that 

=  0  , 

9=0 


(3.68) 


(3.69) 


in  order  to  insure  a  finite  solution  at  the  crown  point. 

It  can  be  observed  that  the  boundary  conditions  for  the  lower 
edge  are  somewhere  between  the  guided-pinned  and  the  guided-c lamped 
cases.  The  degree  to  which  the  boundary  conditions  approach  the 
clamped  case  depends  upon  the  effect  of  the  glued  edges.  For  now, 
only  the  guided-pinned  case  will  be  considered.  From  Figure  3.6,  the 
boundary  conditions  on  the  lowf -  edge  are 


figure  3.6 


CROSS- SECTION  VIEW 
OF  THE  GUIDED- PINNED 

boundary  condition 
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bending  moment  at  (0  =  0)  =  0  ,  (3.70) 

and 

vertical  displacement  at  (0  =  0)  =  0  .  (3.71) 

The  boundary  condition  given  by  (3.70)  will  be  discussed  later 
when  expressions  for  moments  are  written.  The  condition  given  by 
(3.71)  can  now  be  considered.  Thus,  from  Figure  3.7,  one  can  rewrite 
equation  (3.71)  as 

"Is4"“l84'=an5.  (3.72) 

To  consider  the  guided-clamped  case,  one  need  only  replace  equation 
(3.70)  with  the  requirement  that  the  rate  of  change  with  respect  to 
0,  at  the  lower  edge  of  the  shell,  be  zero. 

3.4.3  Finite  Difference  Approximations  for  Derivatives 

The  finite  difference  approximations  to  be  used  in  the  construc¬ 
tion  of  the  model  will  now  be  considered.  For  the  crown  point,  the 
following  difference  expressions  will  be  used: 


J  ^'’j+l  ■ 


(3.73) 


^  ('’j.i  -  +  qj.i> 

and 

where  j  =  0  for  the  crown  point  and  ^0  is  a  half-angle  increment  so 


that 
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0-  =  l/A®  . 


(3.76) 


For  any  interior  point,  the  difference  approximations  used  are 


"j  '  T  ‘"j.!  -  “']> 


(3.77) 


and 


"j  ■  T  <").!  -  *  "j.l'  • 


(3.78) 


For  the  lower  edge,  the  difference  approximations  used  are 


“i  “  f  ‘’j  - 


(3.79) 


and 


"i'  ■  ^  <"1.1  - 

where  j  =  N  for  the  lower  edge.  The  segment  division  and  numbering 
can  be  seen  more  clearly  in  Figure  3.8. 


3.4.4  Stress  Resultants 

The  stress  resultants  in  the  shell  can  be  determined  from  the 
displacements.  This  can  be  done  because  of  the  fact  that  the 
displacements  in  the  shell  are  assumed  known.  From  the  displacements 
one  can  determine  the  strains.  Then,  by  use  of  Hook's  law  for  an 
isotropic  media,  one  can  determine  stresses  as  a  function  of  displace¬ 
ment.  The  stress  resultants  are  then  found  by  integrating  the  stresses 
across  the  thickness.  McDonald  [12],  who  references  Vlasov  [24], 
gives  the  following  equations  for  Che  stress  resultants; 


figure  3.8 


CROSS- SECTION  VIEW 
SHOWING  SEGMENT 
DIVISIONS  AND  NODE 
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e+w^cote)] 


(3.81a) 


8+w  cot  e)] 

(p  80  (pep  0 


(3.81b) 


N-  =  K[u  +w+a(ucot  8  +  v  esc  6+w)] 
6  6  CP 


(3.81c) 


N  =K[a(u  +w)+(ucot6+v  csc0+w)] 
(p  6  (p 


(3.81d) 


where 


2  ' 

12(l-<r^) 


and  K  = 


12(1-<t") 


(3.82) 


The  positive  direction  of  the  stress  resultants  is  shewn  in  Figure  3.9. 
For  rotational  syonetric  vibrations^  one  can  reduce  equations  (3.81) 
to  the  following; 


(3.83a) 


%  -  Dfcrwgg+wgcote] 


(3.83b) 


N  -  K[u  +  w  +  e(ucot6  +  w)] 
6  6 


(3.83c) 


N  =  Krc(u  +  w)  +  u  cot  0  wl  . 
CP  6 


(3.83d) 


The  stress  resultants  of  the  crown  point  can  be  written  frr*ra 
equations  (3.83)  while  noting  the  boundary  conditions  given  by 
equations  (3.67)  and  (3.68).  Hence,  one  obtains 


(3.8i.a) 
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^  K[|  (u^--y  +  'l+cr)w^]  (3.86a) 

and 

=  Kf^  (^1-0  •  O^^eb) 

Cp  J-O  ^  i  o  o 

Now  the  boundary  condition  given  by  equation  (3.69)  can  be  used  in 
connection  with  equation  (3.75)  to  obtain 

w  1  =  3w  -  3w,  +  w_  .  (3.87) 

-1  o  1  2 

Therefore^  inserting  equation  (3,87)  in  (3  85)  results  in 


<"e>J-0  4  f”o  -  *  ”2> 


<3. 88a) 


(Vi=o  “  f“o  -  *  “2I  • 


(3.88b) 


For  an  interior  point,  the  stress  resultants  on  the  k  segment 
are  found  to  be 


<“e>j=k  ■  ? 


+  ci'(Q'  +  2a  cot  — )v,  .  1 

O'  k+1 


D  2  2^ 

(M  )  .  ,  =  7-  [aor  w,  ,  -  2a((jot  -!■  cot  — )w, 
'  cp"j=k  4  k-1  o''  1< 


+  orCao  +  2cot— 


‘%>j=k  -iKZcrcot^- 


+  (l+a)w 


(3.89a) 


(3.89b) 


(3.89c) 


<Vj.k  “  I  *  ‘^“"k+l 


+  (l+a)w,  ]  . 


(3.89d) 


For  the  lower  edge,  i.e.  9  =  0.  the  stress  resultants  are  given  by 


(Me),^N  "  f  [a(Qf- 2acot  e)Wj^_^  -  20(0- a  cot  9) 


*  “""nh  ' 


(3.90a) 


p'd«»v'>luVvv.»';,.. 
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^  [a(cycr- 2cot  e)w^_^  -  ZcyCaor- cot  e)Wj^ 

+  (3.90b) 

(^e)j=N  ^  I  ’*’  (<^‘*"2o-cot  e)'J^  +  2(l+a)Wj^]  (3.90c) 

and 

(Ncj,^j=N  "  f  ^“'^‘^N-1  (o'0^  +  2  cot  e)uj^  +  2(1+(t)w^]  .  (3.90d) 

Hence^  the  boundary  condition  given  by  equation  (3.70)  can  be  used  in 
connection  with  equation  (3.90a)  to  obtain 

"nh-1  "  I  (0?- acot  e)Wj^  -  1  (O'- 2acot  e)w^_^  .  (3.91) 

Therefore^  using  equations  (3.91)  and  (3.72),  where  (3.72)  can  now  be 
written  as 

=  %  tan  0  (3.92) 

one  can  obtain  che  stress  resultant  on  the  lower  edge. 

3.4.5  Equations  of  Free  Vibration 

Lagrange’s  equation  may  be  written  as 

».93) 

where  L  =  T  -  V  is  the  Lagrangian  energy  equal  to  the  difference 
betv/een  the  kinetic  and  potential  energies,  and  where  (')  denotes 
differentiation  with  respect  to  time,  so  that  q^  represents  the 
velocity  of  a  generalized  displacement.  Hence,  for  a  stationery 
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system^  T  =  0  and  V  =  V(q^j  q^^)  so  that  equation  (3.93) 

becomes 


=  0  ,  (3.94a) 


where  q  in  general  represents  u^  v,  and  w  or  3(N  +  1)  displacements 
for  the  discrete  system  being  considered.  For  rotational  symmetric 
considerations,  at  most  2(N  +  1)  -  1  displacements  remain. 

The  total  potential  energy  of  the  discrete  shell  can  be 
obtained  by  summing  the  energies  on  each  segment,  as  now  given  by 
equation  (3.63),  so  that 


N 

V=  +  (up  *  ]  . 

J=0  J  J  J 


(3.94  b) 


The  stretching  and  bending  energies  as  given  by  equations  (3.64)  and 
(3.65)  must  now  be  obtained  across  each  segment  with  appropriate 
consideration  being  given  to  the  boundary  conditions.  Thus,  for  the 
boundary  conditions  to  be  considered  here,  one  can  obtain  2N  simul¬ 
taneous  linear  equations.  In  matrix  notation  this  can  be  written  as 


1-cr 


[S](q)  =  (p) 


(3.95) 


where  [S]  is  the  2N  x  2N  stiffness  matrix, 

[q]  is  the  column  vector  of  displacements, 
ip)  is  the  column  vector  of  external  forces, 
and  h^  has  been  assumed  constant.  The  equations  represented  by 
(3.95)  are  for  the  spherical  cap  with  a  rotationally  symmetric  static 
load.  To  obtain  the  equations  of  free  vibration  for  an  undamped 


The  stiffness  and  mass  matrices  denoted  in  equation  (3.99)  are  net 
presented  in  detail  here.  As  shown  in  Reference  [15],  they  are  found 
in  a  direct  manner  by  integration  of  trigonometric  quantities  to 
obtain  the  equations  of  motion  as  presented  in  Appendix  7.2.  These 
equations  of  motion  are  then  written  in  matrix  form  resulting  finally 
in  equation  (3. 99) . 
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3.5  Acoustic  Radiation 

The  black  box  technique  is  used  to  solve  the  acoustic  radiation 
portion  of  the  problem.  A  computer  program  developed  by  Hess  [8] 
is  used.  For  any  arbitrary  three-dimensional  body,  one  must  input  a 
surface  geometry,  surface  normal  velocity,  and  wave  number  to  be 
considered.  The  input  geometry  is  one  that  represents  a  system  of 
quadrilaterals  to  approximate  the  surface  of  interest.  The  normal 
surface  velocities  are  then  the  velocities  of  each  quadrilateral  making 
up  the  surface.  The  output  is  then  made  up  in  part  by  the  surface 
pressures  corresponding  to  each  of  the  quadrilaterals  and  the  far 
field  pressures. 

3.6  Model  Construction 

3.6.1  General 

As  discussed  previously,  the  model  consists  of  three  basic 
parts.  They  are  the  disk,  the  shell,  and  the  acoustic  radiation 
problem.  These  three  problems  have  now  been  discussed  in  sections 
3,3,  3.4  and  3.5  respectively.  The  manner  in  which  these  problems 
are  now  connected  depends  upon  what  assumptions  one  chooses  to  make. 

The  model  construction  will  now  proceed  in  the  following  order:  shell, 
acoustic  radiation  loads,  and  disk.  The  assumptions  made  will  be 
discussed  at  the  point  at  which  they  are  made  and  are  basically  the 
same  as  those  used  by  Royster  [21], 
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3.6.2  Shell  Vibrations 

3. 6. 2.1  Model  Analysis.  On.  can  recall  that  the  equations  of 
motion  for  an  undamped  spherical  cap  undergoing  free  vibration  were 
given  by  equation  (3.99)  as 

([S]  -  X^[M])(q*}  =  0  .  (3.99) 

This  represents  the  spherical  cap  alone  with  guided  pinned  lower  edge. 
The  physical  transducer,  on  the  other  hand,  does  not  let  the  shell 
vibrate  freely.  The  piezoelectric  disk  provides  resistance  to  the 
free  motion  of  the  shell.  Hence  to  solve  the  true  problem,  one  must 
couple  the  disk  to  the  shell  before  the  free  vibration  problem  is 
solved.  However,  if  the  disk  adds  little  stiffness  to  the  shell  then 
one  can  assume  that  the  solution  of  equations  (3.99)  for  the  eigen¬ 
vector  form  without  including  the  effects  of  the  disk  will  be 
approximately  equal  to  the  solution  with  the  disk  consideration. 

Thus,  making  this  assumption,  one  can  solve  the  eigenvalue  problem 
given  by  equation  (3.99)  by  the  techniques  outlined  by  both  Royster 
[20]  and  McDonald  [13]. 

To  obtain  the  equations  of  motion  for  the  forced  vibration 
problem,  one  considers  equation  (3.95)  where  both  external  forces  and 
inertia  forces  are  to  be  considered.  Thus,  instead  of  equations 
(3.96),  one  considers 

X  =  -phgU  +  X  (3. 101a) 

and 

Z  =  -ph  w  +  Z 
^  s 


(3. 101b) 
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2TTjihg  N 


T ^  [S]  E  b  (tp  )  +2Trpa^h  [M]  Z  b  [cp  )  =  2naJ[p]  .  (3.105) 

1— c  ^rr  ss  -Tr  s 

r=l  r=l 

T 

Premultiplying  by  [cp^}  and  use  of  orthogonality  as  can  be  noted  from 


Reference  [6]  gives 


2TT|ih 


2TTagicpr)'^(p) 


(cpj.)'^CS](cp^)  +  2npagh^b^(cp^]’^[M](cp^) 


Now  from  the  free  vibration  problem,  one  has 


[s]  =  x;[M]  , 


so  that,  noting  equation  (3.100),  one  can  obtain  from  (3.106) 


(3.106) 


(3. 107) 


V  X  2,  1 

b  +  uj  b  =  -r— 

r  r  r  ^ 


(3. 108) 


s  (cpj.)  [M](cp^) 


The  steady  state  solution  for  (3.108)  is  easily  obtained  so  that  the 
steady  state  expression  for  (3.104)  is  given  by 


j(i)t  N 
(q)  =  -  Z 


1 


(9J  ip" 


P^s  r=l  2,  \2. 

^  [!-(—)  ] 


i'Pj.) 


(3. 109) 


where 


*  jiut 
p  =  p^e-" 


and  ^l-M]  (cPj.) 


(3.  110) 


3.  6.  2. 2  Mobility  and  Impedance.  .4s  stated  in  section  3.  6.  2.1, 
the  shell  has  been  assumed  to  vibr-j?*’  independent  of  the  disk.  Hence, 
the  shell  problem  is  solved  .re  . oupled  to  the  disk  problem, 

r 
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considered  in  section  3.3,  through  the  arbitrary  specific  impedance 
applied  to  the  edge  of  the  disk.  Thus,  it  now  becomes  necessary  to 
obtain  the  impedance  around  the  lower  edge  of  the  shell. 

Mobility  is  defined  as 


mobility 


velocity 
force  ’ 


(3.111) 


so  that  from  Figure  3.10  one  can  obtain  the  mobility  at  the  lower 
edge  to  be 

u 

1 ^  .  (3.112) 

F(2na^  sin  0  cos  6) 


Likewise,  one  can  obtain  the  specific  mobility  to  be  given  by 


u„  a  sin  0  [sin  0  sin  l/ot  -  (1  -  cos  l/ o')  cos  01 
Ns 

F  cos  0 


(3. 113) 


To  obtain  the  velocity  in  the  tangential  direction,  u,,  needed  for 

N' 

(3.112)  and  (3.113),  one  can  obtain  from  equation  (3.109)  that 


jujt  N 


iq)  =  -j 


I 


r.l  ,\2, 
U,[l-(— )  ] 


icp^/ 


(3.114) 


The  N  componen'.:  of  equation  (3.114)  then  gives  the  needed  velocity. 
Hence,  the  mobility  or  specific  mobility  for  the  edge  of  the  shell  can 
now  be  obtained.  It  should  be  noted  that,  with  a  knowledge  of  the 
free  vibration  problem  and  the  forcing  iuaction,  one  need  only  con¬ 
sider  the  first  couple  of  terms  of  equation  (3.114)  to  obtain  a  value 
of  mobility  around  the  lirst  resonant  frequency. 
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Impedance  is  the  inverse  of  mobility.  Therefore,  from  equations 
(3.112)  and  (3.113),  one  can  determine  that 


and 

zj  =  1/mJ  .  (3.116) 

3.6.3  Acoustic  Radiation 

For  a  given  shell  geometry,  normal  velocities,  and  wave  number, 
Hess'  program  can  be  used  to  determine  both  the  surface  and  far-field 
pressures.  It  should  be  noted  that  these  values  are  per  unit  density 
of  the  fluid  media.  Two  media  will  be  considered  in  order  to  compare 
with  the  experimental  results. 

The  first  medium  considered  is  air.  Due  to  the  density  of  the 
air,  in  comparison  to  that  of  the  shell,  one  can  neglect  the  resulting 
acoustic  loads.  Thus,  the  force  vector  In  equation  (3.103)  reduces 
to  only  one  component.  Hence,  the  impedance  of  the  shell  can  be  found 
by  using  equations  (3.112),  (3. IIA)  and  (3.115;. 

The  second  medium  to  be  considered  is  water.  The  forced  vector 
now  becomes  complex.  Assuming  that  the  mode  shape,  including  an 
external  load,  is  not  significantly  changed  from  that  of  a  vacuum, 
enables  one  to  greatly  simplify  the  problem.  Tlic  reactive  part  of  the 
acoustic  impedance  can  now  bo  converted  to  an  equivalent  mass.  This 
mass  must  then  be  added  to  the  mass  of  the  original  structure,  and  the 
eigenvalue  problem  resolved  for  the  necessary  data  required  in  Hess' 
numerical  program.  This  iterative  process  can  then  be  followed  until 


the  system  reaches  aqui librium.  Then,  use  of  the  conservation  of 
energy  in  connection  with  the  energy  dissipated  enables  one  to  obtain 
the  final  resistive  part  of  impedance  to  be  applied  to  the  disk. 

3.6.4  Piezoelectric  Disk 

As  discussed  in  section  3.6.3^  with  or  without  significant 
acoustic  loads  one  can  now  obtain  the  resistive  and  reactive  compo¬ 
nents  of  specific  impedance,  i.e.  Q  and  F,  to  be  applied  to  the  disk. 
Then,  varying  the  excitation  frequency  of  the  disk,  one  can  determine 
the  electrical  impedance  of  the  particular  transducer  under  considera¬ 


tion. 
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4.  RESULTS 

The  mathematical  model  of  the  Class  V  f lextensional  underwater 
acoustic  transducer  as  described  in  section  3.6  was  programmed  for 
evaluation  on  the  IBM  360-75  computer  [16].  Results  for  a  particular 
geometry  of  the  transducer  can  be  found  by  defining  the  necessary 
physical  parameters  of  the  transducer. 

For  this  investigation,  experimental  and  mathematical  model 
results  are  compared  for  a  transducer  similar  to  that  shown  in 
Figure  1.2.  The  thin  piezoelectric  disk  has  a  diameter  of  1.5  inches, 
a  thickness  of  0.04  inches,  and  the  material  is  PZT-4  (PZT-4  is  a 
Clevite  Coiporation  trademark  for  a  specific  type  of  piezoelectric 
material).  The  aluminum  spherical  cap  has  a  1,563  inch  inside  radius 
and  is  0. 1  inches  thick.  The  edge  thickness  of  the  spherical  cap  for 
each  experimental  transducer  varies  and  is  noted  with  its  data.  On 
the  other  hand,  the  mathematical  model  is-  restricted  to  constant 

thickness  caps  so  that  no  variation  in  edge  thickness  is  considered. 

2 

The  cxperiraortal  rer-  its  '  j^re^ented  in  Table  4.1  ana  in 
Figures  4.1  and  4.2.  The  mathematical  model  results  in  air  are  pre¬ 
sented  in  Figures  4.3  and  4.4.  Figure  4.3  plots  the  magnitude  of  the 
electric  impedance  in  ohms  versus  the  driving  frequency  in  kilohertz 
(khz)  found  using  the  mathematical  model.  Figure  4.4  shows  both  the 
analytic  and  experimental  results  for  air  in  terms  of  decibels  (db) 
versus  driving  frequency.  The  experimental  results  shown  in  Figure  4.4 
are  for  transducer  number  2. 


^Ibid. 
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Table  4. 1  Experimental  results 

Common  geometric  parameters  of  the  transducer  are  listed  below: 

Piezoelectric  disk:  1.5  inch  diameter 

0.04  inch  thickness 
PZT-4  material 

Spherical  caps:  1.563  inch  inside  radius 

0. 1  inch  thickness 

45°  taper  at  edge  to  edge  thickness 


Transducer 

Number 

Edge 

Thickness 

(inch) 

Frequencies  (khz) 

Air 

Water 

mmamm 

Anti-resonance 

1 

.015 

21.8 

22.6 

12.5 

2 


.030 


22.8 


23.7 


VERSUS  DRIVING  FREQUENCY 


800  ka 


I  2  3  456789  10 

fdO^hz) 

FIGURE  4.3  ELECTRICAL  IMPEDANCE  VERSUS 
DRIVING  FREQUENCY  FOR 
ANALYTIC  MODEL  IN  AIR 


FIGURE  4.4  VOLTAGE  ON  TRANSDUCER  IN  DB  VERSUS  DRIVING 

FREQUENCY  FOR  ANALYTIC  AND  EXPERIMENTAL  RESULTS 
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In  Figure  4.3,  the  analytic  results  of  the  free  disk,  the 
transducer  model  with  a  guided. pinned  boundary  condition,  and  the 
transducer  model  with  a  guided-clamped  boundary  condition  are  pre¬ 
sented.  These  results  show  the  disk  to  be  dominated  by  the  shell,  as 
was  assumed  in  section  3. 6. 2.1.  The  fundamental  resonant  frequencies 
for  the  guided-pinned  and  guided-clamped  boundary  conditions  respec¬ 
tively  are  app.-oximately  20.5  khz  and  33.08  khz,  while  the  anti¬ 
resonant  frequencies  are  approximately  21.5  khz  and  33.1  khz.  Compar¬ 
ing  these  frequency  bounds  to  the  frequency  results  listed  in  Table 
4, 1,  it  is  seen  that  the  experimental  frequency  results  do  lie  between 
the  frequency  bound.®  supplied  by  the  mathematical  model.  It  can  also 
be  seen  that  the  guided-pinned  model  provides  results  much  closer  to 
the  experimental  results. 

In  Figure  4.4,  the  mathematical  results  in  air  for  the  two 
bounding  boundary  conditions  are  compared  to  the  results  of  the 
experimental  transducer  number  2  for  the  electrical  circuit  shown.  As 
noted  above,  the  fundamental  frequency  of  the  experimental  transducer 
is  bounded  by  the  mathematical  results.  However,  away  from  resonance, 
the  experimental  results  are  not  bounded  by  the  results  given  by  the 
two  boundary  co.iditions  in  the  analytic  model.  This  difference  can  be 
attributed  to  several  possible  causes.  First,  the  inaccuracies 
inheren':  in  the  mathematical  model  itself.  These  inaccuracies  can 
result  from  the  uncoupling  of  the  shc’l  and  disk,  from  the  finite 
difference  model  used  to  approximate  the  shell,  from  assuming  a 
constant  thickness,  and  from  the  neglecting  of  the  bond  joints  and 
their  dissipation  of  energy.  A  second  cause  ‘s  the  experimental 
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difficulty  in  obtaining  measurements.  Instrumentation  Inaccuracies 
are  always  present,  and  a  variation  of  +  Idb  is  usually  acceptable. 
Another  point  of  Interest  is  the  dependence  of  the  results  on  the 
series  resistor,  as  shown  by  Royster  [19];  therefore,  a  variation  in 
the  resistor,  due  to  manufacturing  tolerance,  would  produce  a  shift  in 
the  results. 

To  consider  the  transducer  in  water,  acoustic  loads  can  be 
obtained  from  the  numerical  program  developed  by  Hess.  The  loads  must 
be  used  in  an  iterative  manner  with  the  main  program  until  the  system 
converges  to  the  fundamental  resonant  frequency.  In  this  case,  four 
iterations  were  used.  This  required  a  total  of  25  minutes  of  computer 
time.  The  resonant  frequency  of  the  analytic  model  was  11.65  khz,  as 
compared  to  12.5  khz  of  the  experimental  transducer  number  1  in  Table 
4.1.  This  analytic  case  used  the  guided-pinned  boundary  condition, 
since  the  results  in  air  showed  the  condition  to  give  results  nearer 
those  found  experimentally.  A  point  of  interest  is  found  in  that  the 
shaved  edge  of  the  shell  and  the  bond  joint  have  opposing  effects. 
While  the  degree  to  which  each  affects  the  results  is  not  known,  it 
might  be  possible  for  them  to  cancel  each  other  out.  This  could  be  a 
prime  factor  in  the  good  agreement  betwee.i  the  guided-pinned  analytic 
results  and  the  experimental  results. 

In  Figure  4.5,  the  magnitude  of  the  electrical  impedance  of  the 
transducer  in  water,  with  the  guided-pinned  boundary  condition,  is 
plotted  against  the  uriving  frequency,  while  in  Figure  4.6,  the  re¬ 
sistive  part  of  the  electric  impedance  for  this  case  is  plotted 
against  the  driving  frequency,  since  insufficient  experimental 


2  3  4  5  6  7  8  9  to  20  30 


t ( khz  ) 


FIGURE  4.5  ELECTkIC  IMPEDANCE  OF  THE 
TRANSDUCER  IN  WATER  VERSUS 
DR’VING  FREQUENCY 


FIGURE  4.6  RESISTIVE  PART  OF  THE 
ELECTRICAL  IMPEDAHCE  IN 
WATER  VERSUS  THE  DRIVING 

FREQUENCY 


results  are  available  to  which  these  analytic  results  can  be  compared, 
a  complete  evaluation  of  this  part  of  the  model  cannot  be  carried  out 
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5.  SUMMARY  AND  CONCLUSIONS 

The  fundamental  resonant  and  anti-resonant  frequencies  in  air  of 
the  experimental  transducers  were  shown  to  be  bound  by  the  guided- 
pinned  and  gulded-clamped  results,  with  the  guided-pinned  condition 
giving  results  within  approximately  5,9  percent  of  the  experimental 
results.  Therefore,  using  only  the  guided-pinned  boundary  condition 
in  water,  a  fundamental  resonant  frequency  of  the  mathematical  model, 
which  was  only  6.8  percent  lower  than  that  found  experimentally,  was 
found.  From  these  results,  one  can  conclude  that  the  mathematical 
model  with  the  guided-pinned  condition  can  be  used  to  approximate  the 
fundamental  resonant  frequencies. 

In  air,  the  electrical  impedance  of  the  transducer,  found  using 
the  mathematical  model  developed  hetain,  can  be  used  in  the  electrical 
experimental  circuit  to  read  the  voltage  in  db’s  across  the  transducer. 
These  mathematical  results  are  about  2  db's  low  for  the  guided-pinned 
case  and  1  db  lov7  for  the  guidec-clamped  cate  as  compared  to  the 
experimental  results  of  transducer  number  2. 

In  water,  the  electrical  impedance  and  resistive  portion  of  the 
electrical  impedance  can  be  obtained  from  the  developed  model,  used  in 
connection  with  the  numerical  program  wf  ;iess.,  for  the  acoustic  loads. 
However,  due  to  the  lack  of  experimental  data  with  respect  to  the 
electrical  impedance,  no  statement  can  be  made  about  the  reliability 
of  this  portion  of  the  model. 

In  conclusion,  it  can  be  said  that  the  mTthe.nati''al  model 
developed  herein  gives  results  that  compare  well  with  exis  .ing  experi¬ 
mental  data.  However,  additional  experimental  data  must  be  obtained 


with  respect  to  water  loads  before  the  validity  of  the  analytic 
results  for  the  loads  presented  herein  can  be  verified. 
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7.  APPENDICES 


7.1  Obtaining  Cylindrical  Equations  of  State 
from  Cartesian  Equations  of  State 

As  given  by  equations  (3.4)  and  (3.5),  the  Cartesian  tensor 

equations  of  state  are  given  by 


S.  .  =  S^..  ,T,  +  d  .  .E  , 

ij  ijK-t  v-L  mij  m 


(3.4) 


D  =  d  ,  .T,  +  e  E 

n  nK.t  nra  m 


(3.5) 


The  tensor  transformations,  as  given  by  Hawkins  [7]  for  contra- 
variant  tensors  of  rank  one  and  two  respectively,  are  expressed  as 


E  =  ^  E'  , 
m  ax'  n  ’ 
n 


(7.1) 


Sx  dx 

m  _  ^  V  rp  I 

k-t  6x’  st  ’ 

s  t 


(7.2) 


where  E  and  T,  .  represent  any  contravariant  tensors:  x  refers  to  the 

old  coordinate  system,  and  x^  refers  to  the  new  coordinate  system  such 

that  X  =  X  (x ') . 
m  m  n 

Now,  in  order  to  transform  equation  (3.4)  into  curvilinear  co¬ 
ordinates,  one  rewrites  equation  (7.2)  as 


^ '  Sx ' 

c'  —  ^  c 

St  ~  dc ,  ij  ^ 


(7.3) 


SO  that  the  curvilinear  strains  are  in  terms  of  the  Cartesian  strains. 
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Thus,  using  equation  (3.4)  in  equation  (7.3)  gives 


Sx '  3x '  „  bx'  bx'. 

Q I  _  _ s  _ t  _E  _  _ s  _ t  ,  „ 

°st  ScT  c)x^  ijk-l,  kJi  ^x7  ScT  mij  m  ^ 


(7.4) 


which  represents  the  curvilinear  strains  in  terms  of  the  Cartesian 
stress  and  electric  field.  Hence,  using  equations  (7.1)  and  (7.2)  in 
equation  (7.4),  one  obtains  the  curvilinear  equation  of  state 
corresponding  to  (3.4)  to  be 


bx'  bx'  „  bx,  bx,  bx'  bx'  bx 

ct-St  E  k.t.  ,  St.  m  , 

^st  -  ST  ^ijU  ^  ^mn  ST  ST  '‘ijm  ^  n  ’ 
ij  mn  ij-'n 


(7.5) 


Similarly,  the  curvilinear  equation  corresponding  to  (3.5)  is 


^ '  bx,  bx,  bx'  -  bx 

n  c)x  mk-f,  bx '  bx'  st  bx  ®nm  Sc  “s 

m  s  t  ms 


(7.6) 


Hence,  equations  (7.5)  and  (7.6)  give  the  curvilinear  tensor  equations 
of  state  for  a  piezoelectric  material. 

One  should  note  that,  in  general,  the  tensor  components  and 
physical  components  are  the  same  only  in  a  rectangular  Cartesian 
system.  Therefore,  as  shown  in  Fung  [5],  the  tensor  and  physical 
components  of  a  contravariant  tensor  of  rank  two  are  related  by 


(no  sum  on  i  or  j) 


(7.  7) 


where  6  and  e  represent  the  physical  and  tensor  components  respective¬ 
ly  and 


bx^ 

^km  dx,'  bx ' 
k  m 


(7.8) 
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boundary  conditions.  Therefore,  these  equations  must  be  derived  for 

each  segment  and  are,  in  general,  different  from  the  general  interior 

k  /,  J, 

equations.  Defining  the  stiffness  coefficient  ^  ,  where  ^  refers 
to  the  displacement  component  by  which  the  coefficient  is  multiplied, 
and  ^  designates  the  particular  equation  from  which  the  coefficient 
comes,  one  can  obtain  the  following  equations  of  motion; 


j  0,  equation  of  motion  with  respect  to  w^  (crown  point) 

^^0  1  ^“l  ^  ^^0  0  ^”0  ^^0  1  ^”1  ^^0  2  ^'*2  ■  ^*^0  0  ^''0  ^  (7.1-) 


j  =  1,  equation  of  motion  with  respect  to  u, 


Utl.  UU  UW.  QW 

(Sii)ui  .  (S^j_)u2  +  (S^o)Wo  +  (S^,)w^ 


(m“  J  )ui  =  0  ,  (7. 13) 


equation  of  motion  with  respect  to  w, 


J.  /oWU.  ,  ,oWW.  ^  ,«ww.  .  ,_ww., 

*  <^2>“2  *  <*10>"0  *  (®12>“2 


W  W  V  w 

+  (S^  3  )W2  -  (Mj  1  ®  > 


(7.14) 


j  =  2,  equation  of  motion  with  respect  to  u^  (or  k  -  2) 


(S“;)u^ .  (S;;)U2  ms;^)u3  .  (s;;;)wi .  (S“-)W2 


u  u 

-  (M2  2)^2 "  ®  ’ 


(7.  15) 


equation  of  motion  with  respect  to  v. 


W  V  w  w  w  w 

"  <"2  3>-3  "  <"2  4>''4  -  ^*^2  '^"2  °  ^ 


(7.  16) 


J 


j  =  k,  k  :=  3,4,. .  .,N-4,N-3 


equations  of  motion  with  inspect  to 


^  ^  ^  <Sk“k-i>Vi 


. -U  w  ,  ...u  u  . 

<\k'”k  -  <“kk>“k  =  °  > 


(7.17) 


equations  of  motion  with  respect  to  w^^ 

...wu,  ^„w  u  .  ^..ww  .  ,„ww  . 

(\kK  ^  <^kk+lK+l  <\k-2^V2  ^  ^\k-PVl 

WU  VW  WLf 

^^kk^^k  ^^k+l^^k+l  *  ^^kk+2^''k  +  2  '  ^\k^\ 


j  =  N-2,  equation  of  motion  with  respect  to  u^^  2(^  =  N-2) 


(7.18) 


^^N.2  N-3  ^“N-3  ^^N-2  N-2^“n-2  ”*■  ^^N-2  N-1  ^ “N-1  *  ^S.2  N-3  ^''n-3 


^^N-2  N-2  ^''N-2  ■  ^^-2  N-2  ^^-2  ~  ^ 


(7.19) 


equation  of  motion  with  respect  to  w 


^V2N-2)V2  ^  ^V2N-l>Vl  ^  ^VZN^^N  ^  ^  V2  N-4  >  V4 
'^^^N-2  N-3  ^''N-3  ^^N-2  N-2  ^''N-2  ^  ^^N-2  N-1  ^''N- 1 


^^-2  N-2  ■ ‘'N-2  ° 


(7.20) 


j  =  N-1,  equation  of  motion  with  respect  to  u  , 

fl-  I 


^^N-1  N-2  ^^N-2  ^  ^^N-1  N-l  ^“N-I  ^  ^^N-1N^“n  ^  ^^N-1  N-2  ^^-2 


^  ^^N-1  N-1  ^''N-1  ■  1  N- 1  ^  *^N- 1  '  ° 


(7. 21) 
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k+1  ~  ^/oi  +  2'iv  cos  ^sin  l/cy  ’ 

^kk-l""  ^^(l+v)  sln--^i^~^^  sin  1/a  , 


(7.29) 


(7.30) 


u  w  2k  2lc 

"<k  =  -2Q'(1+v)  sin— sin  l/a  +  4(l+v)  cos^sin  1/a  .  (7.31) 

For  the  ^  and  equations,  one  has 

Cl  Ll  =  “'t  sl„  6  sl„  1/,. 

-  (1  -  COS  1/a)  cos  9]  -  4avcos  i^^lHsin  1/a 

+  2[-2  sini^^lilsin  1/a  +  ln  tan^^  -  Intan^^]  ,  (7.32) 

a  2a  2a  ^  ^ 


^  2  .  2(N-1)  .  ,  /  „  2(’N-1)  , 

VlN  ""  sin-i— ^sin  1/a  +  2av  cos  :ilil-i4-sin  l/a 

Or  a 


+  [^  +  a(l+'  tane][(l-cos  1/a)  cos  9-  sin9sinl/a] 
-  av[  cos  6  sin  l/a  +  (1  -  cos  l/a)  sin  9]  , 


(7.33) 


^N-iL2  ^  sin--^y~^^  sin  l/a  , 


(7.34) 


^N-1  N-1  ~  -2a(l+v)  sin— — — sin  l/a +  4(l+v)  cos— — sin  l/a  , 

U; 


(7.35) 
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li  2 

=  2or(l+v)  sin^-^^^sin  l/a+~  (^)^[^tan0  [(a 

-  2 

_  V  cot  6)  (a  -  2v  cot  9)  -  a(2o/  -  3v  cot  0)  +  Of  ] 

~  “  2  3“ 

[  (1  -  cos  l/oO  cos  0  _  sin  0  sin  l/a]  +  a  v  [  cos  0  sin  l/a 

2  _  _  _ 

+  (1  -  cos  1/a)  sin  0]  _  ^tan  e[(l  -  cos  l/a)  cos  0  -  sin  0  sin  l/a 

+  In  tani- In  tan^l^]  -^tan0sin^-^^=^sin  l/a 

_  2HJi_tan0cos^-^^^^sinl/a-^tan0[-2sin^-^^^^sin  l/a 
2  cy  2  O' 

+  In  tan - In  tan —s — ^] }  .  (7.3o) 

a  2a 

The  coefficients  of  the  and  equations  are  given  as  the 

following; 

1  h  4 

Sq  Q  =  4(l+v)  (1  -  cos  l/a)  +  Yj  (^)  [  -g-  [  (1  -  cos  l/a) 

+  2  sin  2/asin  l/a] }  ,  (7.39) 

h  „  4 

®0l"  ■  IJ  ^  X  ^  2  sin  2/asin  l/a] 

s 

3  2 

+  — ^  cos  2/asin  l/a  }  ,  (7.40) 

1  h  „  4 

Sq2=^(^)^  (  cos  l/a)  +  2  sin  2/a  sin  l/a] 

s 

3  2 

+  ^  y  cos  2/asin  l/a  }  , 


(7.41) 
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1  h  -  4 

1  /  Sv  2  r  a 


=  8(l+v)  sin  2/asin  1/a +Yj(— )  ^  {  -^[2(1-  cos  l/a) 


3  2 

+  (4  sin  2/a+  sin  4/a)  sin  l/a]  +2a  v  cos  2/asin  l/a 

2  , 

+  ^  [-2  sin  2/asin  l/a  +  In  tan  3/2a -  In  tan-s—  ] )  ,  (7.42) 


,  h  -  4 

^l2  "  “  IT  ^a“^  ^  ^  [(1-  cos  l/a)  +2(  sin  2/a  +  sin  4/a)  sin  l/a] 

s 

3^  2  ^2 

+  ^ —  [3  cos  2/a  +  cos  4/a]  sin  l/a  +  —  [-2  sin  2/asin  l/a 


+  In  tan  3/ 2a  -  In  tajci  1/ 2a] }  , 


(7.43) 


^  h  4  3  2 

^13  ~  TT  ^  ^sin4/asin  l/a+^-^cos  4/asin  l/a)  ,  (7.44) 

s 

,  h  2  4 

$22  =  2(l+v)  sin4/asin  1/a  +  jj  (j^)  (  ^  [  (1  -  cos  l/a)  +  2(  sin  2/a 

“s 

3  2 

+4  sin  4/a  +  sin  6/a)  sin  l/a]  +  a  v*[  cos  2/a  +  2  cos  4/ a]  sin  l/a 

a^  5  1 

+  -y-  [-2(  sin  2/a+  sin4/a)  sin  l/a  +  ln  tan^ —  In  tan^]}  , 


(7.45) 


,  h  2  4 

^2  3  ^  "  IT  ^  X  ^  sin  6/ a]  sin  l/a 

s 

a3^2  1 

+  — = —  [3  cos  Uja+  cos  6/a]  sin  l/a  +  -=-  [-2  sin  4/Q-sin  l/a 


5  3 

+  In  tan^  -  In  tan^  ]  }  , 


(7.46) 


1  h  -  4  3  2 

ww  1  .s.  2ra  .  ,/  ^av  ,/T  .  ,/ 
^2^  "^12  ^  -^sin6/a  +  — ^  cos  6/a]  sin  l/a  . 


(7.47) 


The  general  expressions  where  k  =  3,  4,  ...,  N-4,  N-3  are 


1 J 
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-ww  0/1J.N  ■  2k  .  ,/  1  /^s.2  r  r  2(k-l) 

\k  "  sin  — sinl/or+y^  (— )  I  -j-  [  sin 


a 


12  'a 


a 


+  sin  ^  ]  sin  l/a  +  a^  [  cos  +  2  cos—]  sin  l/a 

a  '  01  a 


+  ^  [-2(  sin  +  sin^)  sin  l/a  +  In  tan 


a 


1  .  2k-3  ,, 

In  tan ])  . 


(7.48) 


1  h  4 
1  ^  s.  2  f  a 


3  2 


2k  ^g,„2(k+l),„,, 

a 

2k  2(k+l), 


=  .  4(_i)M^[sinfl+sin±ii^]sinl/a 
k  k+1  12  ^a  2  a  a  ■'  ' 

s 


+  [3  cos  ^  +  cos  ^ ]  ®^n  l/a  +  ^  [-2  sin  ^ sin  l/a 


^  ^  2k+l  ,  ^  2k-l,. 

+  In  tan -s- - In  tan  ^r-  J  i  > 


la 


2a 


(7.49) 


„ww  1  ,V2  p  a^  .  2(k+l)  ^a^v^  2(k+l),  .  ,/ 

S,  ,  „  =  (— )  [  4-sin— i - ^  +  ^4; — cos  — i - -\  sin  1/a 

kk+2  12  ^a  ^  4  a  2  a  ■'  ' 

s 


(7.50) 


For  u„  „  and  u„  , ,  one  has 
N-2  N-1 


„«  "  O,,,  ^  .  2(N-2)  .  ,/  J_  1  ,  S.2  ,  a  ..  2(N-3) 

V2N.2'  slt.-i5^sl„l/o+^  (— )  (  -j-  [sin__ 

s 


,  .  2(N-2)  ^  .  2(N-1)t  .  3  2r  2(N-3) 

+  4  sin  — - —  +  sin  — ^ - ]  sin  1/a  +  a  v  [  cos - 

a  a  '  a 

+  2  cos  sin  1/a  +  ^  [-2(  sin  +sin  sin  l/a 

a  2  '  a  a 


^  1  .  2N-3  .  ,  2N-7  ,  , 

+  In  tan  -■--  -  -  In  tan  — :rr—  J  J  , 


2a 


2a 


(7.51) 


X 
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-  2o'(q'-  2v  cot  6)  ]  [  sin  0  sin  l/a  -  (1  -  cos  l/a)  cos  6  ] 

+  cos  +  2  cos  sin  l/a-  cot  S[  cos  I  sin  l/dr 

2 

+  (1-  cos  l/or)  sin  9]  +  ^  [-2  sin  iiO^l^sin  l/a  +  In  tan  Hll 

t  a  '  2a 

-  In  tan^l^]  [_2  sini^|::ilsin  1/a  +  In  tan^l^ 

,  ^  2N-3t  ,  a^  ^ 

-  In  tan —j—]  +—  [  (i_  cos  l/a  )  cos  0  -sin  9  sin  1/a 

+  In  tan  j  -  In  tan  )  . 


(7.53) 


7.3  List  of  Symbols 
radius  of  piezoelectric  disk 
radius  of  shell 

scalar  modal  participation  factor 
electric  displacement  vector 
piezoelectric  strain  constant 
electric  field  intensity  vector 
force  per  unit  circumferential  length 
body  forces  per  unit  volume 
thickness  of  piezoelectric  disk 
thickness  of  shell 

amplitude  of  steady  state  electrical  current 
electrical  current 

Bessel  function  of  first  kind  of  order  N 
planar  coupling  coefficient 
Lagrangian  energy 

mobility  at  lower  edge  of  shell  (j  =  N) 
specific  mobility  at  lower  edge  of  shell  (j  =  N) 
moment  per  unit  length  acting  in  6-direction 
moment  per  unit  length  acting  in  cp-direction 
mid-plane  force  per  unit  length  acting  in  0-direction 
mid-plane  force  per  unit  length  acting  in  tp-direction 
column  vector  of  external  forces 

amplitude  of  steady  state  external  forcing  function 
electrical  charge 

amplitude  of  steady  state  electrical  charge 


q 

{q} 

\ 

[S] 

.u 

sE 


t 


ext 


u,v,w 


u. 

1 

V 


P 

X,Y,Z 

x,z 

Y 

..E 


general  displacement 

amplitude  of  general  displacements  for  steady  state  motion 

column  vector  of  displacement 

normalizing  factor 

stiffness  matrix 

strain  tensor 

stiffness  coefficient  in  stiffness  matrix  [S] 
elastic  compliance  matrix  at  constant  electric  field 
kinetic  energy 
stress  tensor 
time 

strain  energy  due  to  stretching 

strain  energy  due  to  bending 

amplitude  of  steady  state  displacement  of  disk 

strain  energy  due  to  external  forces  per  unit  area 

displacement  components  of  spherical  cap  with  respect  to 

r,  9,  0  directions  respectively 

amplitude  of  displacements  for  steady  state  motion 

displacement  in  the  i^^  direction 

potential  energy 

steady  state  voltage  amplitude 

potential  energy  per  unit  volume 

external  forces  per  unit  area  acting  on  spherical  cap 
external  forces  per  unit  area  other  than  inertia  forces 
electrical  admittance 

Young's  Modulus  for  a  constant  electric  field 
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Y^(x)  =  Bessel  function  of  second  kind  of  order  N 

Z  =  electrical  impedance 

Z'  =  arbitrary  specific  impedance  on  boundary  of  disk 

Z^  =  impedance  at  lower  edge  of  shell  (J  =  N) 

Z*  =  specific  impedance  at  lower  edge  of  shell  (J  =  N) 


a 

r 

AG 


T 

c 

nm 


p 


o* 

n 


0) 


=  inverse  of  half-angle  increment 

=  imaginary  part  of  specific  impedance  on  boundary  of  disk 
=  half-angle  increment 
=  Kronecker  delta 

=  permittivity  matrix  at  constant  stress 
=  absolute  temperature 
=  opening  angle  of  spherical  cap 
=  r^^  eigenvalue 

=  shear  modulus  of  shell  material 
=  density 
=  Poisson's  ratio 
=  entropy  per  unit  volume 

=  eigenvector  corresponding  to  eigenvalue 
=  real  part  of  specific  impedance  on  boundary  of  disk 
=  angular  velocity 
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^e  purpose  of  this  investigation  is  to  develop  a  matheDatical  model  for  the 

ffClasS  V  flextensional  underwater  acoustic  transducer. 

The  transducer  is  approximated  through  the  consideration  of  three  distinct  pro- 

bletts.  The  problem  of  a  thin  piezoelectric  disk  with  an  arbitrary  impedance  <» 

its  edge  is  solved  in  tensr  of  Bessel  functions.  The  shell  vibration  problem 

is  solved  using  a  finite  difference  model  to  approximate  the  shell.  The  acoustic 

radiation  problem  is  solved  by  obtaining  the  source  density  distribution  for  a 

system  of  quadrilaterals  representing  the  transducer.  With  the  source  density  of 

each  quadri lateral ,  the  near  and  far  field  pressxures  and  velocities  can  be  found. 

Utilizing' these  three  components,  a  model  is  then  constructed  for  the  transducer. 

A  comparison  of  the  results  from  the  mathematical  model  with  those  obtained  from 

experiaiaits  is  made  in  order, to  validate  the  model. 
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